Coulomb Effects in Tunneling through a Quantum Dot Stack 
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Tunneling through two vertically coupled quantum dots is studied by means of a Pauli master 
equation model. The observation of double peaks in the current-voltage characteristic in a recent 
experiment is analyzed in terms of the tunnel coupling between the quantum dots and the coupling to 
the contacts. Different regimes for the emitter chemical potential indicating different peak scenarios 
in the tunneling current are discussed in detail. We show by comparison with a density matrix 
approach that the interplay of coherent and incoherent effects in the stationary current can be fully 
described by this approach. 
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I. INTRODUCTION 

Resonant tunneling through coupled quantum dots 
(QDs) is an active field of theoretical and experimental 
research (for a recent review see Various nano-scale 
geometries for double QD systems were under investiga- 
tion: lateral gate-defined devices (e.g. in Ref. ^]) or ver- 
tical mesa structures as proposed in Ref. Q. With the 
recent advances in fabrication of stacked self-organized 
QDs where only few stacks can be selected by electronic 
transport measurements an ideal system for studying "ar- 
tificial molecules" was developed HHIl. The current- 
voltage characteristic of such a double QD stack can 
show peaks with very high peak-to-valley ratio J3| due 
to the sharp energetic resonance condition of levels in 
different QDs. In Ref. 0| the authors also report on 
the frequent observation of double peaks in the tunnel- 
ing current which they relate to the effect of Coulomb 
interaction. In the present paper, we give a detailed the- 
oretical examination of this experimental feature and we 
analyze the mechanisms of the current peak broadening 
and height. We use a Pauli master equation model where 
the tunneling to the contacts is treated by tunneling rates 
and the coupling between the QDs is described by Fermi's 
golden rule. The result of this treatment for the station- 
ary current through noninteracting QDs is in agreement 
with recent density matrix calculations. We will show 
that the occurrence of double peaks in the current de- 
pends on the emitter chemical potential, on the tunneling 
rates, and on the Coulomb interaction strength of elec- 
trons inside the QD as well as in different QDs. Different 
scenarios corresponding to different parameter regimes 
will be considered. 

The outline of the paper is as follows: in Sec. |n]we in- 
troduce the model based on a Pauli master equation; 
in Sec. IIIII we present the analytical and numerical re- 
sults for the stationary current through the double QD 



stack, followed by a discussion. In particular, this sec- 
tion is subdivided into the case of noninteracting elec- 
trons (Sec. IIII A|l where we discuss the effect of cur- 
rent peak broadening and height, and into the case of 
Coulomb interacting electrons (Sec. IIII B|l where double 
current peaks under certain conditions occur; conclusions 
are drawn in Sec. IIVI 



II. MODEL 
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Figure 1: Scheme of the quantum dot stack connected to emit- 
ter/collector contact with tunneling rates Te/Cj respectively, 
and tunneling rate Finter between the quantum dots. 



In Fig. n the double QD system is sketched: QDl is 
connected to an emitter contact with a tunneling rate 
Te, QD2 is connected to a collector contact with a tun- 
neling rate Tc, and both QDs are tunnel coupled mu- 
tually with a rate Finter (for the definition see below). 
We consider one spin degenerate single-particle state per 
QD. The many-particle (Fock-) states v of the QD stack 
are characterized by occupation numbers n['^^ G {0,1}, 
where i = 1,2 labels the QD and s =t, I labels the spin 
degree of freedom. Thus v — {n^^^\ , n'^^K/"'''^'' 
16 different Fock states. 



"-2i') gives 
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The emitter/collector contact is assumed to be in local 
equilibrium characterized by Fermi functions n^^*^ (E) = 
1/(1 + exp((£' — lJtE/c)/kBT)) with chemical potentials 
He/c^ respectively, and temperature T. An applied bias 
voltage V = {^E — A'c)/e drives the QD system out of 
equilibrium. Throughout the calculations we choose the 
sign of the elementary charge e > so that V > Q means 
Mb > and / > describes a net particle current from 
emitter to collector. The single-particle energies Ei in 
the QDs are assumed to shift linearly with respect to the 
bias voltage V: EiiV) — Ei — erjiV with the leverage 
factor rji for the i-th QD. The total energy of the state 
is 



Here f2 is the tunneling matrix element and F = F^ + Tc ■ 
The maximum tunneling rate is Fjnter = 4|r2p/F and 
L{x,w) = (1 + (2x/w)'^)~'^ is a Lorentzian function in x 
with a FWHM of w. For convenience we measure rates 
and energies in the same units here, i.e. setting h = 1. 
Throughout this paper we investigate the steady state 
current so that the stationary occupation probabilities 
evaluated by M • P° = are used. In this case the 
current can be evaluated by counting the tunneling events 
over the collector barrier. Summing the corresponding 
transitions probabilities we find 



I = J2 (1 - 5_ 



{N2{v)~N^{v'))P^, (5) 



(1) 



i=l 



n}t' + n)\' is the particle number of the 



where Ni{v) = n"^^ + nf^ 
z-th QD for state v and C^^O is the charging energy for 
state V which is given by 0, 0, S Q 



U{v) 



\Y.[c-%N^i-) (2) 



with the inverse of the capacitance matrix C, where the 
last term subtracts the self-interaction. 
The evolution of the system is described by the Pauli 
master equation 



—P = K-P 



(3) 



for the occupation probabilities P_ = ({Pj^})-^. The ma- 
trix M_ contains the transition rates M^^i between the 
Fock states v' v. Confining attention to single par- 
ticle transitions all elements M^^i are zero if the total 
particle number f^(v) = 'Ylii^iW) differs by more than 
one from N{v'). For a transition which increases the total 
particle number by unity, i.e. an electron enters the QD 
stack via the emitter/collector barrier, respectively, the 



transition rate reads M^, 



E/C 



E/C'l-F 



(AE^jyr) where 



AE^^> = E{v) — E{u') is the difference of total energies. 
For the inverse tunneling processes the transition rate is 
M^iji = F^/c'(l — n^^'^ {AE^i^)). Transitions which do 
not change the total particle number describe tunneling 
between QDl and QD2: Taking into account the broad- 
ening of the transition due to the finite lifetime in the 
quantum dots, Fermi's golden rule provides us with 



1 



F/2 



^{AE^.'Y 

rinter L{AEi,^i ,T) 



(F/2)2 



(4) 



where the additional factors ensure that the tunneling 
process occurs at the collector barrier and the correct 
sign is taken into account for the current direction. 



III. DISCUSSION AND RESULTS 

A. Noninteracting quantum dots: Current peak 
width and height 

If the energies of single-particle states in two different 
QDs are resonant, a current flows through the QD stack 
and a peak arises in the current-voltage characteristic. 
In this section we give an analytical derivation of the 
current for noninteracting QDs, i.e. setting U{v) = in 
Eq. (Q. Then both spin directions decouple and we can 
restrict ourselves to a single spin direction here. For the 
single-particle levels in QDl and QD2 it is assumed that 
El — A*E ^ ksT and jic — E2 ^ ksT such that rip ^1 

and «0. With P" = (^("o,o)'^(lo)'^(o,i)'^(li))^ 
matrix M simply reads 



M = 
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Tc 
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(6) 



with Z = Fintcri(Ai?(i.o)(o,i) I r). The stationary master 
equation can be rewritten in a rate equation for the occu- 
pation numbers ni and n2 which results in the following 
stationary solution 



ni 
n2 



F, 



FbFc + TZ 
The current |SJl then becomes 



Tc + Z 
Z 



(7) 



I = e 



1 1 



n -1 



(8) 



which can be physically interpreted quite easily: the 
three barriers in the QD stack act as a series connec- 
tion of resistors corresponding to the inverse of tunneling 
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rates. After a straightforward transformation of © one 
obtains 



r^rc 



(9) 



with r„ 



1 

Te ' Fc ' Tint. 

of a hnear dependence of the single-particle levels on the 
bias voltage it follows for the energy difference A£' = 
Ei—er]iV—{E2 — eri2V) = e{ri2—rii){V—VR) with the bias 
voltage for the resonance eVr? = E2z:Mx_ Substituting this 
in Eq. © leads to 



Due to the assumption 



I{y)^eT^L \V-Vr, 




(10) 



As expected, the current-voltage characteristic shows a 
Lorentzian peak at the bias voltage Vr. Interestingly, its 
broadening is not just given by the tunneling coupling 
to the contacts F, but by a line width modified by the 

correction factor 



1 



rerc 



This factor provides a sig- 
nificant increase of the FWHM for the current peak if the 
tunneling matrix element \ is larger than the geometric 
mean ^/VeTc of the contact tunneling rates. 
The expression © for the stationary current was also 
found in |11L [T2. 13] where a density matrix approach was 
used. In the Appendix A we show that both approaches 
give identical results in the stationary state if electron- 
electron interaction is neglected. 

To discuss Eq. ^ let us assume that Yq > > ^e- 
the current is limited by the lowest rate rintcr-^(A£', F) 
for bias voltages V far away from Vr. If the bias ap- 
proaches Vr this rate is increased until it equals the emit- 
ter coupling Te- Now, the emitter barrier starts to limit 
the current, i.e. the current flattens close to the reso- 
nance so that the FWHM value of the current is reached 
for larger distances from the resonance voltage Vr. The 
smaller the ratio of r£;/|r2| the broader the current peak. 
The peak (maximum) current is given by /(Vr) = eF^. 
It has a non-monotonic dependence on the collector tun- 
neling rate Tc'- for weak collector coupling the peak cur- 
rent increases linearly, it shows a maximum at Tc = 2\n\ 
and drops as l/Fp for strong collector coupling. (The 
same holds for the emitter coupling.) On the basis of 
the coherent description this was interpreted in Ref. 
as a competition between enhanced tunneling for weak 
collector coupling and destructive interference in the op- 
posite limit. Hence, it is worth to note that one obtains 
the same effect in a Fermi's golden rule treatment where 
no interference is taken into account for. In this pic- 
ture increasing Fc enhances the transport over the col- 
lector barrier but limits the transport between the dots 
by broadening the transition. 



B. Coulomb interacting quantum dots 

Now, our considerations will be extended to the realistic 
situation of Coulomb interacting electrons. We describe 
the interaction of electrons inside QDl and QD2 with 
charging energies e^[C"^]ii = Ui and e^[C~^]22 = C^2, 
respectively. The Coulomb interaction strength of elec- 
trons in different QDs is given by e^[C~^]i2 = J/inter- 
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Figure 2: Linear bias voltage dependence of the single-particle 
energy levels in QDl (771 =0.26, Ei =79.5 meV) and in QD2 
(??2 =0.68, E2 =118.7 meV). Dashed lines: addition energies 
for the second electron (t/ = I7i = U2 =8meV, Uinta ~0)- 



In Fig.Elthe bias voltage dependence of the QD energies 
arc depicted. The zero-bias single-particle energies Ei 
and leverage factors rji were obtained from the calculation 
of the transmission through a two-dimensional geometry 
as depicted in Fig.Qwithout scattering processes Q. In 
order to describe the experiment, we use the geometry 
in Ref. 6] which gives 771 =0.26, Ei =79.5 meV, and 
7/2 =0.68, E2 =118.7 meV. 

Due to the different slopes, the single-particle levels of 
QDl and QD2 (full lines in Fig. intersect at the point 
Ri . Further resonances can occur by considering the dou- 
ble occupancy of the QDs. Assuming equal charging en- 
ergies for both QDs {U — Ui = U2 = 8meV) parallel lines 
with the additional charging energy U emerge (dashed 
lines in Fig. [SJ. The resonance R2 is due to the inter- 
section of the energies of the doubly occupied QDl and 
the singly occupied QD2. The resonance of the doubly 
occupied states R3 occurs at the same bias voltage as for 
Ri if the charging energies for both QDs are equal. 
A current flow through the QD system is possible at the 
resonance energies E{Ri) and i?(i?2) providing that elec- 
tron states in the emitter contact are occupied at these 
energies. (Typically, these resonances occur at energies 
far above the chemical potential of the collector Fermi, so 
that the electrons can always leave the QD system via the 
collector barrier.) In Fig.|31the dependence of the current 
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Figure 3: Contour plot: current I vs. emitter chemical po- 
tential fiE and bias voltage V for U = Ui — U2 ~ 8meV, 
t/inter = 0, Tfi = Tc = 100/ieV, =10^ieV, T = 4.2 K. 
(Isolines correspond to 0.2, 0.4, 0.8 nA) 



on the emitter chemical potential he and the bias volt- 
age V is shown as a contour plot. Here we assume that 
the tunnel coupling between the dots is rather weak, i.e. 
|0| ^ r^;, Tc, so that the occupation of the dots is deter- 
mined by the respective reservoirs. The conduction band 
edge in the emitter is the zero-point of the energy scale. 
Different scenarios for the stationary current depending 
on He arise: 

(i) He < E{Ri): there is only an exponential small 
current due to thermal activated electrons in the 
emitter. 



(ii) E{Ri) < HE < E{Ri) + U = EiR^): in this case 
the current shows a peak due to the resonance Ri 
of the single-particle levels, the second resonance 
i?2 is energetically available for thermally activated 
electrons from the emitter and gives rise to an addi- 
tional current peak at lower bias voltage which in- 
creases with increasing temperature (this scenario 
was experimentally investigated in Q and was used 
to determine the charging energy U) . 

(iii) E{Ri) + U < He < E{R2): no current peak ap- 
pears. In this regime it is possible to add a second 
electron in QDl. Since the coupling between the 
QDs is weaker than the coupling to the emitter con- 
tact, QDl is mostly occupied with two electrons. 
They leave QDl with the energy Ei + U so that 
they cannot fulfill the resonance condition for Ri . 

(iv) He > E{R2): a current peak due to the resonance 
i?2 occurs with twice the peak height than in the 
case (ii) because two electrons contribute to the 
current here. A current peak for resonance Ri is 
suppressed for the same reasons as in (iii). 



2. Dependence on tunnel coupling 

Now we want to study the behavior with increasing tun- 
nel coupling |r2|. In particular we consider the regime 
(iv) of the last section but now we allow the tunneling 
rates to the contacts to be different {Te — 17/ieV, Tc = 
400 HeY). 
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Figure 4: Current I vs. bias voltage V for varying coupling 
between QDs Te = IT^teV, Fc = 400^eV, T = 4.2K, 
U = Ui=U2^ SmeV, f/intcr = 0, HE = 90meV 



Fig.^shows the bias voltage dependent current for vary- 
ing couplings \ ft\ between the QDs. For low one ob- 
serves a current peak ed V = 75mV corresponding to 
the resonance R2 in Fig. El (compare regime (iv) in the 
last section). According to Eq. IjlOl) the current peak 
broadens with increasing |f2| but one has to consider 
the double occupancy of QDl which modifies the cur- 
rent peak width such that the correction factor becomes 

p^fr^ ■ With increasing \ Q\ a, second current 



1 + 2- 



now 

peak at the bias voltage of Ri appears in Fig. 0] Due to 
the stronger coupling between the QDs single occupancy 
of QDl becomes more likely, since the second electron in 
QDl can tunnel into QD2 off-resonantly because of the 
finite linewidth of the levels. Hence, electrons with en- 
ergy El can leave QDl and contribute in resonance i?i 
which generates a current peak at V{Ri). Hence, the 
possibility of the observation of a double current peak in 
the experiment depends mainly on the tunnel matrix ele- 
ment 51 between states in different QDs. Its sensitivity on 
the spatial separation of the QDs within the stack could 
explain the frequent observation of the double current 
peaks 3. 

Note that for > 260/ieV a small third current peak 
at almost 110 mV appears due to the resonance of the 
energies of the doubly occupied QD2 and single occupied 
QDl (out of the bias range depicted in Fig.|21). 
Underlining the corresponding mechanisms which lead 
to the current peaks in Fig. ^ we show the current- 
voltage characteristics for different charging energies Ui 
with fixed U2 = SmeV in Fig. [S^ (the other parameters 
are the same as for the full curve in Fig. . As one can 
see, with decreasing Ui the low bias current peak shifts to 
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different interaction strengths C/mtcr = 0, 2, 4meV the 
current vs. bias voltage is shown in Fig. |S1 
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Figure 5: Current I vs. bias voltage V for different charging 
energies Ui of QDf; U2 = 8meV; a) Fe = 17^teV, Tc = 
400/ieV; b) Fb = 400/ieV, Fc = 17/ieV; other parameters are 
the same as for the full curve in Fig. ^ 



higher bias voltages and the second current peak remains 
at the same bias voltage position: Reducing Ui leads to 
a parallel shift of the dashed line with smaller slope in 
Fig. 121 so that the intersection point R2 shift to higher 
bias voltages and the resonance Ri leaves unchanged. In 
Fig. 03 the values of F^; and Fc were interchanged so 
that a different physical situation emerges. Here, both 
current peaks shift to higher bias voltages by the same 
amount with decreasing Ui, i.e. the resonance R3 is now 
mainly responsible for the current peak at higher bias 
voltage. Due to the weaker collector coupling the double 
occupancy of both QDs becomes more likely. 
Note, that for equal charging energies in QDl and QD2 
the current-voltage characteristics are invariant with re- 
spect to an interchange of F^; and F^. This is due to 
the electron-hole symmetry in the considered QD sys- 
tem. The necessary condition for this symmetry is that 
the resonances occur at energies deep in the Fermi sea 
on one contact side and unoccupied states on the other 
side. 



3. Coulomb interaction between quantum dots 

The self-organized QDs within one stack grown in 0, 
1^ 1^ have a spatial separation of few nm so that the 
Coulomb interaction between electrons in different QDs 
is a crucial ingredient in a theoretical examination. For 
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Figure 6; Current I vs. bias voltage V for different Coulomb 
interaction strengths between QDl and QD2: f/intcr, other 
parameters are the same as for the full curve in Fig. |1] 



Again a double peak appears whereas the distance of 
the maxima diminishes with increasing C/intcr- If QDl is 
filled with one electron the single-particle level of QD2 
is shifted by C/intcr so that the low-bias resonance R2 
moves to higher bias voltages. The resonance i?i stays at 
the same bias voltage since QD2 filled with one electron 
generates a shift of the single-particle level of QDl by 
the same amount C/intor- This gives for the current peak 
separation 



Ay = 



U — [/inter 



(11) 



Another feature can be seen in the dotted curve in Fig. 
if Al/ becomes smaller, the height of the high-bias peak 
increases. This peak structure is quite similar to that 
observed in H. 



IV. CONCLUSIONS 

The non-equilibrium current-voltage characteristic of a 
system of two QDs coupled in series was investigated by 
a Pauli master equation approach. The different leverage 
factors of the energy levels in the QDs leads to resonances 
at certain bias voltages. If the emitter contact provides 
electrons at the resonances a current through the QD 
system is driven. The height and the width of the arising 
current peaks in the current voltage characteristic was 
examined analytically for the case of vanishing Coulomb 
interaction. Here it was found that the width of the 
resonance increases if the tunnel coupling between the 
QDs surpasses the geometric mean of the couplings to the 
contacts. These results are in full agreement with a fully 
coherent approach. Including the Coulomb interaction a 
double current peak can arise for asymmetric coupling 
to the contacts and a tunnel coupling of the order of the 



6 



Coulomb interaction strengths. The strong sensitivity of 
the couphng parameters to the barrier thicknesses in self- 
organized QD stacks can provide a possible explanation 
of the frequent observation of double current peaks in the 
experiment. 



ally, probability conservation is demanded: p^^ = 1 
and pcb — pIc holds. Eq. I|A1(I can be transformed in 
a matrix equation of the form u ~ A ' IL with u = 

{Paa, Pbb,Pcc, Pdd, Re{pbc), Im(pbc))^ and 
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Appendix A: DERIVATION OF THE MASTER 
EQUATION FROM DENSITY MATRIX 
APPROACH 

We start with the modified Liouville equation for the 
double QD system (Fig.P) derived in Ref. [13 • The fol- 
lowing abbreviations for the Fock states {f) will be used: 
\a) = |0,0), \b) = |1,0), |c) = |0,1), and \d) = |1,1). 
Then the time evolution of the corresponding density- 
matrix elements is given by 



Paa = —^EPaa + ^CPcc 
Pbb = ^EPaa + ^CPdd + i^{Pbc - Pcb) 
Pec = -rpcc - i^iPbc - Pcb) 
Pdd = —^CPdd + ^EPcc 

1 



Pbc = iAEpbc + i^ipbb - Pec) - -j^Pbe 



with A£; 



El and r = Fb + Fc 



(Al) 
Addition- 



A 



-F 



Te 



V 





Ve 










-n 






Tc 











-iF 
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2f7 


-AS 



(A2) 



In order to derive the Pauli master equation for the oc- 
cupation probabilities P^, — p^jy one has to get rid of 
the non-diagonal elements pbc in (|A1|) by truncating the 
6x6 matrix HA2p to the 4x4 matrix ©. This can be 
accomplished by setting the time derivatives of Re(pbc) 
and Im(pbc) to zero. Then one can solve the algebraic 
equations for p^**** = ^(pcc — pbb)/{^E + jF/2) and sub- 
stitute this expression in the equations for pbb and pec- 
This immediately leads to the terms Z — A£2^(r/2 y^ ™ 
© as introduced by Fermi's golden rule in Sec. [J 
The necessary assumption pbe = is justified if one of 
two conditions is met: (i) In the limiting case <C F 
the relaxation of pbe{t) to p^**** occurs on a fast time 
scale. This is the condition for sequential tunneling which 
underlies our derivation of the Pauli master equation, (ii) 
In the stationary case pbc — holds independently of the 
magnitude of Jl. Therefore, the Pauli master equation 
provides reliable results also in the strong coupling limit 
if one restricts oneself to the stationary current. 
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